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METHOD FOR EMBEDDING COLOR IMAGE INFORMATION INTO A BLACK- 
AND-WHITE IMAGE 

FIELD OF THE INVENTION 
The present invention generally relates to methods for converting color 
images to black-and-white images and, more particularly, to methods for 
embedding texture information into a color image's black-and-white reproduction 
to enable the conversion of the black-and-white image back into the color image. 

BACKGROUND OF THE INVENTION 
5 When color images are printed on black-and-white printing devices the 

color information contained in the original image has to be converted or 
otherwise mapped to black-and-white. This has to be achieved in a manner 
wherein sufficient color information comprising the color image is retained in the 
derivative black-and-white image such that the image itself is not lost in the 
10 recomposition. Since regions of different color having a similar lightness look 
about the same, subtle color differences in the original color image can be 
difficult to distinguish in the resultant black-and-white image. In some instances, 
this inability to distinguish color differences can compromise the content of the 
original picture and thereby compromise one's understanding of the resulting 
15 graphic or chart. 

Colors are generally mapped to textures in order to impose some degree 
of differentiability other than the density of toner (luminance). Different colors 
should produce different textures which should be more readily discernible. 
Methods approach this problem by mapping colors to one or more texture 
20 patterns. Alternatively, only primary colors are mapped to texture patterns so 
that the combinations of primary colors end up being mapped to combinations of 
texture patterns. 
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As wavelet theory is important to an understanding of the subject matter 
of the present invention, the following additional background is provided. 

It is well known from Fourier theory that a signal can be expressed as the 
sum of a, possibly infinite, series of sines and cosines. This sum is also referred 
5 to as a Fourier expansion. A disadvantage of a Fourier expansion is that it has 
only frequency resolution and no time resolution meaning that, although all the 
frequencies present in a signal might be determinable, one does not know when 
they are present. To overcome this problem, several solutions have been 
developed which are more or less able to represent a signal in the time and 

10 frequency domain at the same time. 

The idea behind time-frequency joint representations is to cut the signal of 
interest into several parts and then analyze the parts separately. It is clear that 
analyzing a signal this way will give more information about the when and where 
of different frequency components, but it often leads to a fundamental problem 

15 as well: How to cut the signal? Herein lies the problem. Suppose that we want 
to know exactly all the frequency components present at a certain moment in 
time. This short time window is cut using a Dirac pulse (/(t)=1 at t=0 and f{\)=0 
for all other t), then transform it to the frequency domain. 

One problem is that cutting the signal corresponds to a convolution 

20 between the signal and the cutting window. Since convolution in the time 
domain is identical to multiplication in the frequency domain and since the 
Fourier transform of a Dirac pulse contains all possible frequencies then the 
frequency components of the signal begin to smear over the frequency axis. 
This is the opposite of the standard Fourier result as we now have time 

25 resolution but no frequency resolution whatsoever. The underlying principle of 
this problem is due to Heisenberg's uncertainty principle which states, in part, 
that it is impossible to know the exact frequency and the exact time of 
occurrence of this frequency in a signal. As such, a signal can not be 
represented as a point in the time-frequency space. 

30 The wavelet transform or wavelet analysis is one solution to overcome the 

shortcomings of a Fourier transform. In wavelet analysis, the use of a fully 
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scalable modulated window helps resolve the signal-cutting problem because, as 
the window is shifted along the signal, the spectrum is calculated for every 
position. As this process is iterated with a slightly shorter (or longer) window for 
every cycle, a collection of time-frequency representations of the signal, all with 

5 different resolutions, a multi-resolution analysis results. 

In the case of wavelets, one does not speak of time-frequency 
representations but rather time-scale representations, wherein scale is the 
opposite of frequency as used for Fourier transforms. The large scale is 
considered the big picture, while the small scales show the details. In a way, 

10 going from large scale to small scale is like zooming in on the details. 



The wavelet transformation preferred herein is the Continuous Wavelet 
Transform (CWT) as described more formally as: 



where * denotes complex conjugation. This equation shows how a function f{\) is 
decomposed into a set of basis functions, called wavelets. The variables s and 

scale and translation, are the new dimensions after the wavelet transform. 
The wavelet transform (1) is calculated by continuously shifting a continuously 
20 scalable function over a signal and calculating the correlation between the two. 
As such, the wavelets themselves are generated by: 



where s is the scale factor, r is the translation factor and s-1/2 is for energy 
25 normalization across the different scales. 

An admissibility condition is used to analyze and reconstruct a signal 
without loss of information is satisfied by the square integrable functions ^(t), as: 



Brief Wavelet Transformation Theory 
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where ¥(0) is the Fourier transform of ^(t). The admissibility condition implies 
that the Fourier transform of ${\) vanishes at the zero frequency, i.e., that 
wavelets must have a band-pass like spectrum. Thus: 

i2 



= 0 

co=0 



(4) 

A zero at the zero frequency means that the average value of the wavelet 
in the time domain must be zero, and thus oscillatory. In other words, !/^(t) must 
10 be a wave. 



J\|/(/)<& = 0 



(5) 

The CWT maps a one-dimensional signal to a two-dimensional time-scale 
joint representation that is highly redundant. The time-bandwidth product of the 
15 CWT is the square of that of the signal and, for most applications, is inefficient. 
Discrete wavelets overcome this problem because they can only be scaled and 
translated in discrete steps. By modifying the wavelet in (2), we get: 

1 t k% qSq 
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(6) 

20 wherein, j and k are integers and Sb>1 is a fixed dilation step. The translation 
factor o depends on the dilation step. 

When discrete wavelets are used to transform a continuous signal, the 
result is a series of wavelet coefficients referred to as the wavelet series 
decomposition. The effect of discretizing the wavelet is that the time-scale 

25 space is now sampled at discrete intervals. Usually s 0 = 2 is chosen so that the 
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sampling of the frequency axis corresponds to dyadic sampling. For the 
translation factor, usually o = 1 in order to have dyadic sampling of the time axis. 

The following illustrates the localization of the discrete wavelets in the 
time-scale space on a dyadic grid. 
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An important issue a decomposition is reconstruction. Signals can be 
reconstructed from its wavelet series decomposition because the energy of the 
wavelet coefficients lies between two positive bounds, such that: 

2 - 2 



4/r*zi(Av,,)i <B\f 



10 (7) 

where ||f|| 2 is the energy of f(t), A>0, B<~ and A,B are independent of /(t). 

The discrete wavelets can be made orthogonal to their own dilations and 
translations, which means: 



1 if j = m and k-n 
0 otherwise 
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An arbitrary signal can be reconstructed by summing the orthogonal 
wavelet basis functions, weighted by the wavelet transform coefficients: 

/(0 = ZtO",*)Vm(0 

0) 

5 which shows the inverse wavelet transform for discrete wavelets. 

In (4), the wavelet has a band-pass like spectrum. From Fourier theory, 
compression in time is equivalent to stretching the spectrum and shifting it 
upwards: 

i (<A 

F{f(at)} = -F - 
a\ {a J 

10 (10) 

In other words, a time compression of the wavelet by a factor of 2 
stretches the frequency spectrum of the wavelet by a factor of 2 and shifts all 
frequency components up by a factor of 2. Using this, the finite spectrum of a 
signal is covered with the spectra of dilated wavelets in the same way our signal 
15 in the time domain is covered with translated wavelets. To get good coverage of 
the signal spectrum, the stretched wavelet spectra should touch each other, as 
shown. This can be arranged by correctly designing the wavelets. 
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As a result, a series of dilated wavelets can be seen as a band-pass filter 
20 bank. The ratio between the center frequency of a wavelet spectrum and the 
width of this spectrum is the same for all wavelets and is the fidelity factor Q of a 
filter and, in the case of wavelets, is a constant-Q filter bank. 

When the signal has infinite energy it will be impossible to cover its 
frequency spectrum and its time duration with wavelets. Usually this constraint is 
25 formally stated as: 
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This means that the L2-norm of signal f{\) is finite. Thus, natural signals 
normally have finite energy. 
5 If the scaling function is thought of as just a signal with a low-pass 

spectrum, it can be decomposed into wavelet components and expressed as: 

9(0 = ZY(M)M^(') 

(12) 

Since the scaling function <P(t) is selected in such a way that its spectrum 
10 is neatly fitted in the space left open by the wavelets, (12) uses an infinite 
number of wavelets up to a certain scale f. 

scaling function spectrum (tp) 
y 1 c . ork . ...wavelet spectra 

* — | k V 
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This means that if a signal is analyzed using the combination of scaling 

function and wavelets, the scaling function by itself takes care of the spectrum 
15 otherwise covered by all the wavelets up to scale y, while the rest is done by the 

wavelets themselves. Thus, the number of wavelets reduces to a finite number. 
By introducing the scaling function, the problem of having an infinite 

number of wavelets is circumvented and a lower bound for the wavelets is set. 

Of course a scaling function, instead of wavelets, results in information being 
20 lost. However, from a signal representation standpoint, no information was lost 

since it is still be possible to reconstruct the original signal as will be shown. But, 

from a wavelet-analysis standpoint valuable scale information gets discarded. 

The width of the scaling function spectrum is therefore an important parameter in 
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the wavelet transform design. The shorter its spectrum the more wavelet 
coefficients and the more scale information. But there are limitations on the 
number of wavelet coefficients one can handle. However, in the Discrete 
Wavelet Transform (DWT) this problem is more or less automatically solved. 
5 The low-pass spectrum of the scaling function allows the admissibility 

condition to be re-stated, similar to (5), as: 

J(p(/)^/ = l 

(13) 

This relationship shows that the 0-th moment of the scaling function 

10 cannot vanish. (If the degrees of freedom in a wavelet transform design are not 
only used on creating vanishing moments for the wavelet, but are equally 
distributed among the scaling function and the wavelet to create vanishing 
moments for both functions, we have Coiflets.) 

Therefore, if one wavelet can be seen as a band-pass filter and a scaling 

15 function is a low-pass filter, then a series of dilated wavelets together with a 
scaling function can be seen as a filter bank. 

If the wavelet transform can be regarded as a filter bank, then wavelet 
transforming a signal can be regarded as passing the signal through this filter 
bank. Analyzing a signal by passing it through a filter bank is known as sub- 

20 band coding. The outputs of the different filter stages are the wavelet transform 
coefficients and the scaling function transform coefficients. 

The filter bank needed in sub-band coding can be constructed in several 
ways. One is to build many band-pass filters to split the spectrum into frequency 
bands. The advantage of this is that the width of every band can be chosen in 

25 such a way that the spectrum of the signal to analyze is covered in the places 
where it might be interesting. The disadvantage is that every filter has to be 
designed separately. Another is to split the signal spectrum into two equal parts, 
a low-pass and a high-pass part. The high-pass part contains the smallest 
details (granularity) of interested and we stop here. We now have two bands. 

30 However, the low-pass part still contains some details and therefore can split it 



again and again, until we are satisfied with the number of bands created. In this 
way, an iterated filter bank is created. The advantage is that we have to design 
only two filters. The disadvantage is that the signal spectrum coverage is fixed. 
The process of splitting the spectrum is illustrated graphically as: 
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What remains after the repeated spectrum splitting is a series of band- 
pass bands with doubling bandwidth and one low-pass band. (Although in 
theory the first split produced a high-pass band and a low-pass band, in reality 
the high-pass band is a band-pass band due to the limited bandwidth of the 

10 signal.) The same sub-band analysis is performed by feeding the signal into a 
bank of band-pass filters of which each filter has a bandwidth twice as wide as 
his left neighbor (the frequency axis running left to right) and a low-pass filter. 
The wavelets give the band-pass bands with doubling bandwidth and the scaling 
function provides the low-pass band. The wavelet transform becomes a sub- 

15 band coding schema using a constant-Q filter bank (although not restricted to 
constant-Q filter banks). If the wavelet transform is implemented as an iterated 
filter bank, the wavelets do not have to be specified explicitly. 

In order to use the results achieved with a discrete signal, the wavelet 
transform must be made discrete (discrete wavelets are not time-discrete, only 

20 the translation and the scale step are discrete). 

Given, in (12), a scaling function expressed in wavelets from minus infinity 
up to a certain scale /, if a wavelet spectrum is added to the scaling function 
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spectrum, a new scaling function is produced with a spectrum twice as wide as 
the first. The first scaling function can therefore be expressed in terms of the 
second because all the information required is contained in the second scaling 
function, such that: 

<p(2- / /)=5;A J+ i(*)<K2- /+1 /-*) 

5 * 
(14) 

This relationship states that the scaling function, at a certain scale, can be 
expressed in terms of translated scaling functions at the next smaller scale, i.e., 
in more detail. Since the first scaling function replaced a set of wavelets, the 
10 wavelets can be expressed in terms of translated scaling functions at the next 
scale. The wavelet at level / becomes: 

k 

(15) 

Since signal f{\) can be expressed in terms of dilated and translated 
15 wavelets up to a scale f{\) can also be expressed in terms of dilated and 
translated scaling functions at scale /, as: 

/(/) = 2>,(Ar)<p(2^-Ar) 

k 

(16) 

If, in (16), we step up one scale to wavelets have to be added in order 
20 to keep the same level of detail. 

/«) = £ Vi (*)<P(2^' - k) + £ Yj-, (kM2 J -* t - k) 

k k 

(17) 

If the scaling function (pj, k (t) and the wavelets ^j, k (t) are orthonormal then 
the coefficients Aj.i(k) and Yj.^k) are the inner products: 



10 



Vi(*) = (/(/XVj.*(0> 

(18) 

If <P j>k (t) and #j, k (t) are replaced by suitably scaled and translated versions, 
e.g., replace 2't in (14) and (15) by ^t-k, and arrange terms: 

5 

X J _,{k) = Y J Krn-2k)X J {m) 

m 

(19) 

m 

(20) 

10 These state that the wavelet coefficients and scaling function coefficients, 

on a certain scale, can be found by calculating a weighted sum of the scaling 
function coefficients from the previous scale. 

Recall that the scaling function coefficients came from a low-pass filter 
and recall how a filter bank was iterated by repeatedly splitting the low-pass 

15 spectrum into a low-pass and a high-pass part. The filter bank iteration started 
with the signal spectrum. So, if we imagine that the signal spectrum is the output 
of a low-pass filter at the previous (imaginary) scale, then our sampled signal can 
be regarded as the scaling function coefficients from the previous (imaginary) 
scale. In other words, our sampled signal f(k) is equal to A(k) at the largest 

20 scale. 

A discrete weighted sum, as in (19) and (20), is the same as a digital filter. 
Since the coefficients Aj(k) come from the low-pass part of the splitted signal 
spectrum, the weighting factors h{k), in (19), form a low-pass filter. And since 
the coefficients Aj(k) come from the high-pass part of the splitted signal 
25 spectrum, the weighting factors g(k), in (20), form a high-pass filter. This means 
that (19) and (20) together form one stage of an iterated digital filter bank and 
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the coefficients h{k) become the scaling filter and the coefficients g(k) the 
wavelet filter. In (1 9) and (20), the scaling and wavelet filters have a step-size of 
2 in the variable k. The effect of this is a sub-sampling property that only every 
other Aj(k) is used in the convolution with the result that the output data rate is 
equal to the input data rate. The sub-sampling property solves the problem of 
how to choose the width of the scaling function spectrum. 

Every time the filter bank is iterated the number of samples for the next 
stage is halved. Thus, you are left with just one sample (in the extreme case). 
This is where the iteration stops. This also determines the width of the spectrum 
of the scaling function. Normally the iteration stops where the number of 
samples becomes smaller than the length of the scaling filter or the wavelet filter, 
whichever is the longest. So the length of the longest filter determines the width 
of the spectrum of the scaling function. As such, (19) and (20) can be 
implemented as one stage of an iterated filter bank, such as: 





Thus, the highly redundant continuous wavelet transform with its infinite 
number of unspecified wavelets is reduced to a finite stage iterated digital filter 
bank easily implemented on a computer. Redundancy is removed by using 
discrete wavelets. A scaling function solves the problem of the infinite number of 
wavelets needed in the transformation. 

The reader is encouraged to peruse one or all of the following. Kaiser, G., 
A FRIENDLY GUIDE TO WAVELETS . Boston: Birkhauser, 1994, Weiss, L. G., 
WAVELETS AND WIDEBAND CORRELATION PROCESSING , IEEE Signal 
Processing Magazine, Jan 1994, pp. 13-32. Sheng, Y., WAVELET 
TRANSFORM, in: Transforms and Applications Handbook . A.D. Poularikas (ed). 



12 



pp. 747-827, Boca Raton, Fl (USA): CRC Press, 1996 (The Electrical 
Engineering Handbook Series). Burrus, C.S., Gopinath, R.A., Guo, H., 
INTRODUCTION TO WAVELETS AND WAVELET TRANSFORMS, A PRIMER , 
Upper Saddle River, NJ (USA): Prentice Hall, 1998. Hubbard, B. Burke, THE 
5 WORLD ACCORDING TO WAVELETS , Wellesey, MA (USA): A K Peters, 1996. 
Mallat, S.G., A THEORY FOR MULTIRESOLUTION SIGNAL 
DECOMPOSITION: THE WAVELET REPRESENTATION , IEEE Transactions on 
Pattern Analysis and Machine Intelligence, Vol. 11, No. 7 (1989), p. 674-693. 
And, Vetterli M. and C. Herley, WAVELETS AND FILTER BANKS: THEORY 
10 AND DESIGN . IEEE Transactions on Signal Processing, Vol. 40, No. 9 (1992), 
p. 2207-2232. 

BRIEF SUMMARY 

A method is disclosed for converting color images to textured 
monochrome images such that image regions with similar luminance but different 

15 chrominance look different when converted to black-and-white to preserve color 
information therein. This texture-encoded color information can be extracted at a 
later time to reconstitute the original color image. The present method involves 
first converting the color image to a luminance-chrominance colorspace. The 
chrominance data is decomposed into 4 channels of chrominance. A wavelet 

20 transformation of the luminance channel is preformed. Scaled chrominance 
channels are then mapped to a number of wavelet sub-bands. Once mapped, 
the wavelet transform is inverted to generate textures proportional to the original 
colors of the color image having an amplitude proportional to the chroma of the 
original color. The black-and-white image is embedded with this texture 

25 information into the image's grayscale component. 

In order to retrieve the color image from the textured grayscale image, a 
wavelet transform regenerates the scaled chrominance channels mapped 
therein. Those wavelet sub-bands comprising embedded chroma are extracted 
and scaled to the size of the original image. The remaining sub-bands are 

30 inverted to produce the Y-image thereof. Separate images now exist in Y, Cr, 
and Cb. The recombination of the YCbCr images produce the colored image 
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that was previously embedded into the grayscale components of the textured 
black-and-white image. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The preferred embodiments and other aspects of the invention will 
5 become apparent from the following detailed description of the invention when 
read in conjunction with the accompanying drawings which are provided for the 
purpose of describing embodiments of the invention and not for limiting same, in 
which: 

Figure 1 shows a flow diagram of the present method of converting the 
10 color image to a textured grayscale image; 

Figure 2 shows a flow diagram of the present method of reconstructing 
the color image from the textured grayscale image; 

DESCRIPTION OF THE SPECIFICATION 
A method is presented to convert color images to textured monochrome 
15 images such that image regions with similar luminance but different chrominance 
look different when converted to black-and-white in order to preserve color 
information within the black-and-white image. 

With reference now being made to Fig. 1, the present method involves 
first converting an original color image 10 from its colorspace into a luminance- 
20 chrominance colorspace comprising separate images of Y (luminance), Cr 
(chrominance red/green) and Cb (chrominance blue/yellow), shown at 12, 14, 
and 16 respectively. If the original colorspace was, for instance, RGB then the 
following conversion is preferably utilized to map the RGB colorspace into YCbCr 
space: 

25 Y = 0.299R + 0.587G + 0.1 14B 

Cb = - 0.1687R - 0.331 3G + 0.5B + 128. 
Cr = 0.5R-0.4187G-0.0813B + 128. 
After mapping, the images of the Y, Cr, and Cb are each the same size as the 
original. In other words, if the original image 10 is size 400x400 pixels, then 
30 each of the luminance-chrominance images (12, 14, 16) would be size 400x400 
as well. 
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The Y-image 12, undergoes a wavelet transformation thereof to produce 
18. Preferably, a quick transformation is used to map 8-bit images into 8-bit 
wavelets. If an 8-bit representation results in a blurring of lines, then a 16-bit or 
32-bit is preferred. The wavelet transform separates various frequencies 
5 representative of information comprising the original color image. Note that the 
data of 18 is grouped into sub-bands (represented in the figure as blocks) and 
the frequency of the sub-bands generally increase from left-to-right and from top- 
to-bottom. In other words, the data of the upper-left sub-band 18A! represents 
the lowest frequency while the data in 18D represents the highest frequency. 

10 Each of the sub-bands identified as 18AL 18A 2 , 18A 3 , and 18A 4 , in the upper-left 
quadrant of 18 of Fig. 1, comprise the lower frequency data with each sub- 
quadrant thereof representing spatial information of the scaled-down original 
color image. The data contained in sub-band 18B in the upper-right quadrant 
represents vertical edge information. The data contained in sub-band 18C in the 

15 lower-left quadrant of 18 represents horizontal edge information. And, the data 
of sub-band 18D in the lower-right quadrant represents diagonal edge 
information. 

The size of the data array comprising the wavelet transformed luminance 
image is the overall same size as the original image in pixels. If the original color 

20 image is, for instance, 400x400 pixels then the size of the wavelet transformed 
data array is also 400x400 pixels. Each of the sub-bands 18Ai, 18A 2 , 18A 3 , and 
I8A4, in the upper-left quadrant are 1/16 the area of the original image. In this 
instance, each sub-band therein is 100x100 pixels each. Each of the other 
quadrants 18B, 18C, and 18D are 200x200 pixels each. 

25 Reference is now being made to 20 and 22 of Fig. 1, wherein the Cr- 

image 14 of size 400x400 and the Cb-image 16 of size 400x400 are scaled down 
so that each can later be mapped to quadrants of the transformed data of 18. 
Cr-image 14 is thus scaled to 1/4 the area to get Cr 20 comprising a 1/2 scaling 
in the vertical and a 1/2 scaling in the horizontal. The same is done to Cb-image 

30 16 to get Cb 22. Such a scaling can be achieved by taking, for example, every 
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other pixel in the vertical and every other pixel in the horizontal. Other ways of 
scaling an array to 1/4 size are well known in the arts. 

Reference is now being made to 24, 26, 28, and 30 of Fig. 1 . Once the 
scaling has been done to Cr-image 14 and the Cb-image 16 to form scaled 
5 arrays 20 and 22, these values are separated into positive and negative 
components such that: Cr+ = Cr*u(Cr); Cr- = -Cr*u(-Cr); Cb+ = Cb*u(Cb); Cb- = - 
Cb*u(-Cb); where u(x) is a unit step function. The step function as used herein is 
such that positive values remain while negative values are zeroed. Thus, the 
Cr+ 24 comprises the Cr values of pixels with positive values and zeroed for the 

10 pixels with negative values. The Cr- 26 comprises absolute values of Cr for 
pixels with negative values and zeros for pixels with positive values. Similarly, 
the Cb+ 28 comprises the values for pixels with positive values and zeros for 
pixels with negative values. The Cb- 30 comprises absolute values of Cb with 
values of pixels with positive values zeroed. 

15 The next step of the present invention involves mapping each 

chrominance channel Cr+ 24, Cr- 26, Cb+ 28, and Cb- 30 to sub-bands as 
shown. The number of sub-bands replaced with chrominance channels is 
minimally two, one for Cb and one for Cr although a single sub-band can be 
used for a hue correlate (since it is usually hue not chroma which carries the 

20 most information). Alternatively, three are mapped; one each for Cb and Cr and 
one for plane C = Cb*u(Cb) + Cr*u(-Cr). The C plane can be incorporated as an 
extra plane to replace yet another wavelet sub-band. Four sub-band 
replacements are preferred when the desire is to convert back to color at a future 
time. 

25 As described in Fig. 1 , the chrominance channels 24, 26, and 28 do not 

require any further scaling as their size corresponds to the size of the sub-band 
data which they replace. Channel 24 replaces the sub-band at 18 as 18B. 
Channel 26 replaces the sub-band at 18 as 18D. Channel 28 replaces sub-band 
18C. Channel Cb- 30 is to be mapped to a wavelet sub-band which is 1/4 the 

30 area. Again, this reduction can be achieved by a scaling of 1/2 in the vertical 
and 1/2 in the horizontal although other scaling techniques are known in the arts. 
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Thus, Cb- channel 30 is further scaled to 1/4 size, shown at 32, to fit the sub- 
band it replaces. The result of the sub-band replacements is shown at 34. Note 
that 3 of the 4 sub-bands comprising the upper-left quadrant, namely, those 
corresponding to 18A 1( 18A 2 , and 18A 3 retain their original values. Only 18A 4 is 
5 replaced by the data of Cb- at 32. Alternatively, the sub-bands that are replaced 
by the various chrominance channels can be changed. Similarly, any of the 
channels can be rescaled and replaced at 32. Also, wavelets with more levels of 
decomposition can be used. This technique is not restricted to a two-level 
wavelet. 

10 Once the sub-bands have been replaced by the chrominance channels, 

an inverse wavelet transform 36 is applied to generate textured image 38. 

Attention is now directed to Fig. 2. To recover the color image from 
textured image 38, a wavelet transformation 40 thereof produces wavelet sub- 
bands collectively at 41 . This essentially corresponds to 34 of Fig. 1 although a 

15 small error may creep into the data depending on the precision of the method of 
calculation used. The Cb+, Cb-, Cr+, Cr- channels are extracted. Thus, sub- 
band Cr+ (41 B) is extracted to Cr+ at 42. Sub-band Cr- (41 D) is extracted to Cr- 
at 44. Sub-band Cb+ (41 C) is extracted to Cb+ at 46. Sub-band Cb- (41 A4) is 
extracted and scaled to become Cb- at 48. 

20 The absolute values of Cb+ Cb- Cr+ Cr- are taken and used to recompose 

Cb = |Cb+| - |Cb-| and Cr = |Cr+| - |Cr-|. Thus, channel 42 is combined with 
channel 44 to recompose Cr-image 50. Likewise, channels 46 and 48 are 
combined to recompose Cb-image 52. Then each of these channels is scaled to 
match the size of the original image. Thus, Cr 50 is scaled to Cr-image 56 and 

25 Cb 52 is scaled to Cb-image 58. 

With reference now being made to 54, once the sub-bands have been 
extracted, these sub-bands are set to contain all zeroes. The other remaining 
sub-bands corresponding to I8A1, 18A 2 , and 18A 3 still contain data. At 60, an 
inverse transform of this data in 54 produces Y-image 62. 

30 Post-processing is preferably applied to luminance Y and chrominance Cb 

Cr to remove noise, increase saturation, improve contrast, or sharpen, as the 



process of half-toning, printing, and scanning tends to degrade these 
characteristics. 

The recomposition of components Y-image 62, Cr-image 56 and Cb- 
image 58 form the image, shown collectively at 64, now in YCbCr colorspace. 
This colorspace can be subsequently mapped to RGB or another colorspace, if 
desired. 

While particular embodiments have been described, alternatives, 
modifications, variations, improvements, and substantial equivalents that are or 
may be presently unforeseen may arise to applicants or others skilled in the art. 
Accordingly, the appended claims as filed and as they may be amended are 
intended to embrace all such alternatives, modifications variations, 
improvements, and substantial equivalents.. 
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